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The semi-classical limit of quantum-mechanical scattering in two dimensions (2D) is developed. 
We derive the Wentzel-Kramers-Brillouin and Eikonal results for 2D scattering. No backward or 
forward glory scattering is present in 2D. Other phenomena, such as rainbow or orbiting do show 
up. 
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I. INTRODUCTION 

The quantum mechanical scattering problem in two 
space dimensions (2D) has been addressed Refs. 0,0, 0]. 
The inherent circular symmetry in 2D renders the free 
partial-wave radial function to be the ordinary Bessel 
(regular at the origin) and Neumann (irregular) func- 
tions, in contradistinction to the the three-dimensional 
(3D) case where these functions become the spherical 
counterpart (spherical Bessel and Neumann functions). 
2D scattering has application in surface physics, besides 
being an interesting pedagogic topic at the graduate level. 

The semi-classical (SC) limit of the 3D scattering ad- 
dressed in Refs. 0,0 is an important subject as it shows 
the connection between quantum and classical mechanics 
in a particularly transparent form. It has also important 
application in atomic, molecular and nuclear scattering 
where, at large energies, one usually encounters sums 
over prohibitively large numbers of partial waves. Such 
sums are converted into integrals and efficiently evalu- 
ated in SC scattering. Further, the use of semi-classical 
ideas is a common practice in the description of wave op- 
tics phenomena such as the meteorological rainbow and 
glory [§[ . The concepts employed in the wave optics field 
found a very natural adaptation in SC 3D scattering, 
where one deals with the optics of matter waves. 

It would be only natural to extend the SC limit to 2D 
scattering, which should find interesting application in 
surface physics where the results of 2D quantum scatter- 
ing are often applied. Here we develop the SC version of 
2D scattering using the deflection function [j| approach. 
The SC limit covers a situation where the potential varies 
very little over a distance of the order of de Broglie wave 
length. The potential V need not be weak so long as en- 
ergy E ^> | V\ . Hence the domain of validity is different 
from that of quantum Born approximation [7[. Under 
these situations the semiclassical path concept is valid 
and the quantum wave-function is replaced by the semi- 
classical Wentzel-Kramers-Brillouin (WKB) wave func- 
tion 

Here we develop the SC version of 2D scattering using 



the deflection function [j|[8[ approach. We present results 
for 2D scattering in the WKB and Eikonal approxima- 
tions Q as well as a comprehensive discussion of rainbow 
and glory scattering and orbiting. From the asymptotic 
form of 2D SC wave function we extract the phase shift 
and construct the Eikonal result for scattering (cross sec- 
tion) using this phase shift. The same Eikonal result can 
also be obtained in 2D from an evaluation of the quan- 
tum Born approximation to scattering amplitude using 
the SC WKB wave function, as has been done in 3D pj. 
We leave this as an exercise for the reader, which can be 
performed in a straightforward fashion following the 3D 
analysis of Ref . [7| ■ 

In Sec. UH we give a brief account of 2D quantum 
scattering following Ref. 2]. In Sec. Ill II a self-contained 
discussion of SC 2D scattering is presented. Specifically, 
we present results for WKB wave function and phase 
shift, Eikonal phase shift and scattering amplitude, SC 
cross section, and rainbow cross section. In Sec. IIVI a 
summary of our findings is given. In Appendix A we 
present a study of Rainbow scattering cross section at 
different angles where we also present a project for the 
reader. In Appendix B we have left an exercise. 



II. QUANTUM SCATTERING IN 2D 

The time-independent Schrodinger equation for re- 
duced mass fj,, potential V(r), energy E, describing scat- 
tering is 

AvV(r)+^(r) = %). (1) 
In polar coordinates [r = (r,9)], Eq. (Q]) becomes [HQ 



19/ dtp 



r dr \ dr 



1 d 2 jj 



+ [k 2 - U(r)]ib = 0, (2) 



where k 2 = 2^E/h 2 and U{r) = 2fiV(r)/h 2 . Here U(r) 
is a central potential (depending only on the magnitude 
of r). Equation (2) can be solved using the separation 
of variables method [H, 0], namely, one writes ip(r) = 
Y, m R m{r)T m {e) with d 2 T m {0)/d9 2 = -m 2 T m (9), so 
that T m (8) = 7T2 cos(m8), thus yielding 
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Introducing R m = R m ri , Eq. Q can be rewritten as 
d 2 R m (k,r) 



dr 2 



F(r) 



+ F(r)R m (k,r) = 0, 
[m 2 -\ 



U(r)- 



(4) 
(5) 



The function R m (t, r) should be a regular at r = with 
asymptotic form similar to that of the free solution (no 
potential). Setting U(r) equal to zero in Eq. one 
gets the ordinary Bessel equation, whose solutions are the 
regular Bessel [J m (kr)] and irregular Neumann [N m (kr)} 
functions, with asymptotic forms for kr ^> 1 [f|: 
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The asymptotic form of radial function R m (kr), being a 
linear combination of J m (kr) and N m (kr), is taken as |2j 



Rrn(k,r) > A m (fc) 



fcr 



mi: 



+ S m {k) , (8) 



where 6 m (k) is the scattering phase shift and the constant 
A rn is given by A m = 2i m e m (2n)~h iS ™ ( V ,with e m = 2 
for m/0 and eo = 1. 

The scattering amplitude is written in terms of 5 m (k) 
as 
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m=0 



cos(m6»)e m e' i5m(fc) sin5 m (fe). (9) 



The differential scattering cross section is da2D/d9 = 
\f(k,6)\ 2 /k and the total cross section is <J2D{k 2 ) = 
J„ 27r \du2D I d9)d9 = 4 e m sin 2 8 m (k). In 2D the scat- 
tering amplitude has the dimension of L 1 / 2 (L in 3D) and 
cross section has the dimension of L (L 2 in 3D). 



III. SEMI-CLASSICAL SCATTERING IN 2D 

When the local de Broglie wave length is short com- 
pared to the distance over which the potential changes 
appreciably, concepts of classical scattering can be used 
to describe quantum scattering. The problem is formu- 
lated in terms of classical quantities, such as the impact 
parameter, b, and the deflection function, Q(k 2 ,b) 0, Q: 
the validity criterion is 89 opt /Q(k 2 , b) <C 1 for all b and E. 
Here 89 op t is the optimal angular dispersion of the scat- 
tering particle around a classical path. This corresponds 
to the high-energy region where many angular momenta 
contribute and angular momentum is often treated as a 
continuous variable and the quantum-mechanical partial- 
wave sum is replaced by an integral. 



The useful functional relation of the deflection function 
0(fc 2 , b), obtained from the Hamilton- Jacobi equation of 
motion, in both 2D and 3D is given by 0], 



Q(k 2 ,b) = ir-2bk / 

J rn 



drr 2 F 2 ( r ) 



(10) 



where ro is the classical turning point: F{tq) = 0. 

At this stage we recall that Eqs. (JH) and can also be 
solved by the SC WKB approximation 5] with the ansatz 
R = (- J F 1 )- 1 / 4 exp[/ r ' o dr'i-F) 1 / 2 H near r = 0. As we 
cross the turning point and reach the asymptotic region 



we take (-F) 
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In this fashion the wave- function is real in both regions. 
Then for r > r [1 



kr - 



dr'F 1 ' 2 



driF 1 / 2 - fc) - fcr 



(11) 
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From Eqs. |(5J| and (fT2"]) . the WKB phase shift becomes 



(k,m) = rmr/2 + I dr(F 1/2 - fc) - fcr .(13) 



cWKB, 
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Differentiating Eq. (fT3"|) with respect to m = kb we get 



d5Y D KB (k, 



= m 



dm 2 
Comparing with Eq. PH|) we establish 
d5l™ B (k,m) 



drr- 2 F- 1/2 . 
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dm 



(14) 



(15) 



At very high energies many angular momentum con- 
tribute and the so-called Eikonal approximation can be 
derived by replacing the m-sum in Eq. ^ by an integra- 
tion with the identification m = kb [5j: 



j-Eik 
J 2D 



(k,b) 
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(2tt) 1 /2 



db cos(kb9) 



(16) 



The Eikonal phase shift 6^ is now obtained from Eq. 
(fT3|) by expanding F in a power series in U(r)/k 2 , setting 
r = b, m — kb and m 2 — 1/4 w fc 2 6 2 , and keeping the 
lowest order term: [^] 



rEik 
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2fc 
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(17) 



Incidentally, in 3D the Eikonal approximation is similar 
to Eqs. (fTB"|) and p?|) except for the replacement db — > 
{2Tr) 1 ' 2 bdb, cos(fc66») -> J (kb9) [|. 

We now recast the scattering amplitude into a sum of 
integrals [5| using the Poisson sum formula 0: 



K — — OQ 



dyF{y)e 2 ^. (18) 
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Formula (fl~8]) can be applied to scattering amplitude Q 
once cos(to#) is expressed in terms of exponentials: 



/(M)=/ (+) (M) + / l_; (M) 



(-)i 



(19) 



where the amplitudes f^ + \k,8) and f('(k,6) corre- 
spond to the exp(im#) and exp(— imO) branches of the 
cosine function, respectively. Then using Eq. (|18|) 



oo 



x / exp [i{25(k, to) + 2-Kmn ± m9}] dm, (20) 
Jo 

where now the discrete sum over angular momentum is 
replaced by an integration and the SC phase-shift S(k, to) 
is essentially the WKB phase shift. (We drop the label 
WKB in the following.) 

The integrals in Eq. (|20p can be evaluated in the SC 
regime of stationary phase approximation (SPA) when 
the phase shift, being an action divided over Ti, is very 
large. The formula of the SPA is based on evaluating in- 
tegrals of the type I = J dmexp[i£(m)], where the phase 
£(to) = [25(k, m) + 2-7TKTO ± m9] is real on the line of in- 
tegration and analytic in some region surrounding it. If 
£(to) varies rapidly with to, the integral will be vanish- 
ingly small, except in cases where £ has extrema on the 
line of integration. One thus expands ((m) around the 
extremum and keeps terms up to second order in to. For 
one extremum point, m s 



C(to) = C(m s ) + C"(m.)(m - to s ) 2 /2, 



(21) 



where prime denotes derivative. The integral can now be 
performed to yield 

/oo 
exp [i("(m s )(m - m s ) 2 /2] dm 
-OO 

= [2 m /C'(m s )Ve^\ (22) 

where the limits of integration have been extended to 
±oo since the contribution of the far away points are 
strongly suppressed by oscillation. The above result can 



be generalized for several extrema, to 



U). 



I = 



E 
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(23) 



The power of the SPA is that it allows replacing an infi- 
nite sum, such as the to-suhi in / by a sum of only few 
contributions arising from the stationary points. The 
Poisson sum can now be performed, by seeking the sta- 
tionary points of the total phase C( m ) y i a dC t {m)/dm = 
0:, viz, 



±6 = -2 



d5{k, to) 
dm 



2k7T. 



In Eq. (|24p re represents the number of times the particle 
circles around the scatterer for an attractive interaction. 



This phenomenon is known as orbiting. Taking re = 
in Eq. (j2"4"|) and comparing with Eq. ([15]) we find that 
the scattering angle 6 is related to the deflection function 
0:0 = =f6* [5j, where the minus sign refer to and the 
plus sign to f('. If O is positive f(~^ dominates and vice 
versa. The case = corresponds to a situation where 
/( + ) and are comparable in size and the cross section 
should show a Fraunhofer-type interference. From Eqs. 
(HI, (EU), and dHJ) we get 



/(±)(fc } 0) = V" e i{±m< ^ >e+2S+7r/4) 

(±) 



dB 
dm 



(25) 



in 



(±) 



where the TO-derivative is evaluated at the stationary 
points to 

da 2D 



The differential cross section now is 



dB 



= ll\f i+) + f { - ] \ 2 \ 



E 



J-rClassical 
aa 2D 



d6 



(26) 

7 (±) (fc 2 ,0),(27) 



where (da§p ' / d&) (±j is the classical differential 
cross section arising from stationary point to s (±) 



(28) 



" j^Classical " 

aa 2D 
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de 


d6 
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dm 



and the interference terms I^(k 2 , 0) are associated with 
the two branches of the amplitudes. 

If only one stationary point contributes to, say, f(~ 
then the SC differential scattering cross section is just 
the classical one (no interference term). This is similar to 
the 3D Coulomb scattering Q where the exact quantum 
mechanical cross section is identical to the classical one, 
as the deflection function 0(to s ) is a monotonic function 
of the classical impact parameter, identified as b = km s . 
In the 2D case, this identification is just the same and 
one finally obtains for the 2D classical differential cross 
section 



J-rClassical 
aa 2D 

d9 



de 



db 



(29) 



to be compared to the classical 3D differential cross sec- 
tion [| 



J-rClassical 
aa 3D 

dn 



dQ 



db 



(30) 



If the differential cross sections are written as differentials 
with respect to 0, (dil = 2~ n'm6d6.) we have from Eq. 



^—Classical 

^ = 2 *b 

de 



djnb 2 ) 

de 



(31) 



(24) to be compared to the 2D cross section (|2"9")l , namely 



J-rClassical 
aa 2D 

de 



db 



(32) 
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where in Eqs. (J3TJ) and (J32J) we have set |9| = \6\ 
as the stationary point condition requires. Correspond- 
ingly, the total cross sections integrated over all angles 
are agassicai = 7r6^ ax and CT ^ assical = 26 max , where the 
integration is limited to a maximum |&| = 6 max - 

From Eq. ([29]) we see that when the deflection function 
is stationary with respect to variations of impact param- 
eter b, e.g. dQ/db — 0, the differential cross section is 
infinite. Moreover, Eq. (|10p shows that 0(b) will not 
be a monotonic function of b and the differential cross 
section may exhibit interference like in a rainbow ||. 
Here angular deflection, as measured by the deflection 
function, can be positive or negative. The experimental 
scattering angle 9 is limited by < 9 < tt. Hence the val- 
ues of impact parameter b, for which O > 0, contributes 
more to f(~^ and negligibly to and vice versa. If 
dQ/db = for positive O, will show rainbow char- 
acteristics and vice versa. If deflection function 0(to), 
with m = bk, is expanded about the rainbow value m r of 
angular momentum m r = kb r at which Q'(m r ) = 0, 



9(m) = 0(m r ) + 6"(m r )(m - m r ) 2 /2 + 



(33) 



then the phase £ of Eq. (|2TH) becomes £± = 26 (m r ) + 
0(m r )(m — m r ) + 0" (m r )(m — m r ) 3 /6 + 2%Km ± m9, 
[compare with Eq. (|21[) ]. (For repulsive potential there 
is no orbiting and k should be set equal to zero. How- 
ever, the results for amplitude and cross section be- 
low do not depend on n.) Then the resulting integral 
in Eq. ([2"0|) has the same form as in Eq. ([2"2"|) and 
can be expressed in terms of the Airy function 0, Q 
Ai(x) = (1/2tt) exp[i(te + t 3 /3)]dt, leading to 



l/(M)l = 



'2tt 



[Q"(m r )/2]V : 



■Ai 



— Or 



[9"(m r )/2] 1 /3 



and the cross section is 



da 2 p 
d9 



2 5 / 3 fc 1 / 3 ^ 

(|rf 2 e/rf6 2 | br ; 



-Ai 2 



2 1 / 3 k 2 / 3 (e-e r ) 



(\d^e/db 2 \ br 



(34) 



(35) 



where 9 r = 0(m r ) is the rainbow angle. The Airy func- 
tion being an oscillating function the rainbow-like inter- 
ference pattern in the differential cross section (j35|) is 
confirmed. For 9 < 9 r , the expression ([33]) oscillates and 
the particle is scattered into the illuminated region. For 
9 > 6 r , the Airy function dies out rapidly, as the particle 
scatters into the dark side of the rainbow. This feature is 
of a purely quantum nature related to the phenomenon 
of tunneling; classically, there is no scattering for 9 > 6 r , 
and thus the name " dark" . Quantum mechanics supplies 
some "illumination" in this classically forbidden region. 

Besides the rainbow, one has the effect of glory scat- 
tering in 3D, both near 9 = (forward glory) and near 
= tt (backward glory). For example, when back scat- 
tering is possible for impact parameter other than zero 
(b 0), one has 9 = tt and cross section ([30| can be infi- 
nite due to the vanishing of sin 9. It occurs whenever of 
Eq. (|10[) goes smoothly through or through an integral 



multiple of tt [|| . In 2D cross section ([29]) the sin 9 term 
is absent and one cannot have glory scattering. However, 
one can obtain oscillating cross section near 9 = or tt 
owing to near equality of and f(~' [described by 
Bessel function(s)]. 

The phenomenon of orbiting will take place in both 2D 
and 3D when the effective potential U(r) + b 2 k 2 jr 2 has 
a maximum for a certain bo and the energy of the parti- 
cle coincides with this maximum One has to verify 
whether d(U + b 2 k 2 /r 2 )/dr = for a certain bo and ro 
and check if k 2 = U(ro) + k 2 b\jr\. From expression (fl0|) 
for , one can show that in the case of orbiting Q 



e(k 2 ,b) 

e{k 2 ,b) 



const + cln[(6 — b )/b ], 
const + 2cln[(6 — b)/b ] 



b>b (36) 
b < b . (37) 



As b — > bo, Q(k 2 ,b) can be infinitely large and for at- 
tractive interaction the particle can go around the center 
many times. Then the differential cross section made 
from many components — one from each orbit — could 
be greatly different from the classical results ([29)1 and 
(|3H)l and will show interference effect [5j . 

The interference terms of the SC 2D cross section ([27]) 
is absent in cross section ([29|l for a single stationary point 
n s . The case of two stationary points in Eq. ([27]) leads 
to an interference term and deserves special attention. 
We assume that deflection function Q(k 2 ,b) is a double- 
valued function. For a given 9, there are two impact 
parameters, b\, b 2 {n\ = kb\ and n 2 = kb 2 ), that con- 
tribute to Eq. ([25]) . If there are two contributions to , 
say, f(~) in Eq. ([2l>]) . then dQ/dm will be positive at one 
of these points, and negative at the other. This results in 
an extra phase of tt/2 in the contribution to the ampli- 
tude in Eq. (|25[) coming from the latter stationary point. 
Then from Eq. ([23]) the scattering cross section can be 
written as 



d6 
sin[fc(&i 



db x 



de 

K 



de 



db\ db 2 
d9 d9 



(38) 



Equation Q38p clearly shows the wave nature of the scat- 
tering process. The interference term involving the sine 
term is purely quantum mechanical would be absent in 
classical scattering. The classical scattering cross sec- 
tion is an incoherent sum of contributions. The inter- 
ference pattern has an angular period, Pg (angular dis- 
tance between two adjacent maxima in the scattering 
cross section), which can be read off from the above, as 
Pe = 27r/|fc(6i — 62) j . When b\ approaches b 2 , the anal- 
ysis leading to Eq. ([55]) breaks down and one has to 
resort to a more elaborate treatment, called the Uniform 
Approximation [lO], which provides a generalization of 
Eqs. (|34p and ([3l)|) involving the Airy function and its 
derivative. 
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IV. SUMMARY 

We have presented a comprehensive discussion of semi- 
classical (SC) scattering theory in two dimensions. We 
have presented and discussed classical and SC cross sec- 
tion in 2D, WKB and Eikonal approximations. Al- 
though, rainbow scattering and orbiting take place in 2D, 
there is no forward and backward glory scattering. 
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APPENDIX 



A. An Application of Rainbow Scattering 



The semiclassical cross section (135]) exhibits the essen- 
tial quantum effects of tunneling and interference. To see 
this more clearly we use the asymptotic form of the Airy 
function Ai(x) as x — > ±00 [ll| 



Ai(a;) 
Ai(x) 



exp(-2a; 3 / 2 /3) 
2V^ 1/4 : 
sin[2(-cc) 3 / 2 /3 + tt/4 

V^(-^) 1/4 



x > 0, 



(39) 

x < 0. (40) 



Using the asymptotic form (|3"9"|) . we find that the cross 
section (|3"5t exhibits a decay in angle that goes like 



da 



2D 



exp 



4fcV2 (0-8 r f 2 

'-\„r) 1/2 



3 (| ge 



(19 



(41) 



in the classically forbidden angle region, 9 > 9 r , which is 
reminiscent of tunneling effect [classically, there is zero 
scattering in this dark side of the rainbow, see Eq. |29p]. 

In the classically allowed, illuminated, region, 9 < 9 r , 
the cross section exhibits interference arising from the 
contributions of the two stationary phase points. Using 
the asymptotic form (|4*U|) in Eq. (|55|) . the oscillation in 
cross section is given by 



do- 



2^811 



2D 



2fcV2 (e r -ef 2 
3 i^eii/2 



+ 



(id 



1 rf 2 e 1 



(42) 



One then finds the period of these oscillations as follows. 
Call the argument of the sin 2 function above, g{6) = 
^ (^-g) 3/2 + * Suppose ate = 9 1 < 9 r , we have 



(H§Fm i/2 



a maximum in cross section implying sin 2 [<7(#i)] 
and 3(6*1) = 7r/2. At a slightly smaller angle, 9 — 
9±, one has another maximum implying sin 2 [9(6*2) 



= 1, 
h < 
= 1 



corresponding to 5(02) = tt/2 + tt. To obtain an estimate 
of the local period of oscillation Pg of cross section (|42p , 
defined by Pg = 9\ — 62, we expand g(9) near 9 — 9\ as 
g{9) = g(9i) + (9 - 9 1 )g'(9 1 ). Consequently, the period 
Pg at 9 = 9\ is Pg 1 = —ir/g'(9i). Hence the generic local 
and angle-dependent period for any 9 < 9 r is 



Pn 



dg(e) 



kV2 



(I^M 1/2 



(43) 



Thus the smaller the angle 9 compared to the rainbow 
angle 9 r , the smaller will be the period. 

The angle 9 — 9\ is called an Airy maximum; in the 
meteorological rainbow, it is the last maximum in the 
Airy function, just below 9 — 9 r , and corresponds to the 
primary bow, whereas the other maxima are the super- 
numeraries. The second maximum at an angle slightly 
smaller than that of the primary bow is called the sec- 
ondary bow. One usually sees both of these bows in the 
clear sky after rain, the secondary bow being quite dim 
with the order of colors inverted compared to the more 
familiar primary bow. 

The period formula (|43[) is an approximate one valid 
under condition (|33|) used in deriving the rainbow scat- 
tering (|55|) . Keeping only two terms in Eq. (|3"3")l . one 
obtains the solution for two turning points mi and mi as 
m 1{2) =m r ± V2V0r-9/\O"{m r )\ 1 / 2 (recall, 0(m r ) = 
9 r , 9(m) = 9, and 9 r > 9). Thus the period of the Airy 2 
function in Eq. (|35p is Pg = 2n/\mi — m^l, which is the 
same as that derived in Eq. (|4"3"| (note that m = kb). 

The correct period should be obtained by staying at 
angles much smaller than the rainbow angle [maximum 
of O(m)] and using the semiclassical approximation em- 
ployed in obtaining Eq. (|25p and the quadratic approxi- 
mation of the phase £, given by Eq. (f2"Tj) . [In the semi- 
classical approximation leading to Eq. (pj)) . Eqs. (j2"2")l 
and (|23p diverge at the rainbow caustic. The ampli- 
tude (j3~4")l is free from that divergence but the period 
of oscillation (|4"3"| obtained from this amplitude is not 
quite right.] The uniform approximation of Ref. [To| 
combines the good features of both approaches and pro- 
duces the correct period of oscillation. However, the 
algebra is quite involved and without proof we quote 
the result for the angle-independent period so obtained: 
Pg = 1-k l\m\ — TO2I, where mi and m-2 are the m values 
at two stationary points. The period (|43p can be shown 
to be equivalent to the result so obtained after using the 
relations mi = kb\ and m.2 = kbi- (We leave this as a 
project for an advanced reader.) 



B. An Exercise for Reader 



We consider an interaction of the form 



.4 


r > R c , 




r 


A 






Hi 


2R C 





r < R c , 



(44) 
(45) 
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where A and R c are constants. The above potential 
should not be confused with the Coulomb interaction be- 
tween charged particles, which in 2D behaves as ln(r). 
The deflection function for the above interaction shows 
a rainbow maximum and exhibits the phenomenon of 



orbiting. The condition for rainbow is that the energy 
E > 3A/(2R C ). The calculation of the deflection func- 
tion is presented in [![. We leave the calculation of 2D 
SC scattering with this potential as an exercise. 
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